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Compact (ferro- and antiferromagnetic) sigma-models and noncompact (hy- 
perbolic) sigma-models are compared in a lattice formulation in dimensions 
d > 2. While the ferro- and antiferromagnetic models are essentially equiva- 
lent, the qualitative difference to the noncompact models is highlighted. The 
perturbative and the large N expansions are studied in both types of models 
and are argued to be asymptotic expansions on a finite lattice. An exact 
correspondence between the expansion coefficients of the compact and the 
noncompact models is established, for both expansions, valid to all orders on 
a finite lattice. The perturbative one involves flipping the sign of the cou- 
pling and remains valid in the termwise infinite volume limit. The large N 
correspondence concerns the functional dependence on the free propagator 
and holds directly only in finite volume. 
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1. Introduction 



Nonlinear sigma-models with maximally symmetric riemannian target spaces naturally 
come in dual pairs, one compact and the other noncompact. The generalized classical 
spin systems associated with both the compact and the noncompact target spaces have 
a variety of applications, see e.g. [1, 2, 3] for the less familiar noncompact models. 
Since the compact models are much better understood and the target manifolds can 
be related to those of the noncompact models by analytic continuation, it is natural 
to try to relate also the quantum (or statistical mechanics) properties of the compact 
to that of the noncompact models. For definiteness we consider here as target-spaces 
the iV-dimensional sphere and its dual the hyperboloid H N , however the qualitative 
aspects should be the same for other dual pairs. 

The goal of this note is to present two examples of such relations for invariant correlation 
functions. The first one concerns perturbation theory, the other the large N expansion. 
Provided one has chosen a formulation where the expansions are valid asymptotic ex- 
pansions, it turns out that knowledge of the expansion coefficients in the compact model 
allows one to infer those in the noncompact model. However it is crucial that the 
asymptotic expansions are known to exist beforehand in both systems independently. 
The lattice formulation is especially suited to address this and we shall see that on a 
finite lattice the relevant asymptotic expansions do exist; for the perturbative one in 
section 2 below and for the large N expansion in a separate paper [4] . The perturbative 
correspondence simply involves a sign flip of the coefficients and in a formal expansion 
(dimensional regularization and minimal subtraction) has been noted to low orders in 
[5, 6] and in the literature on Riemannian sigma models. Here we present a proof of the 
correspondence to all loop orders, initially on a finite lattice; but the termwise infinite 
volume limit should exist on account of the expected lattice counterpart of David's the- 
orem [7, 8]. The perturbatively defined correlation functions, viewed as functions of the 
lattice points and of 1//3, are thus likewise related simply by flipping the sign of f3. 

The large N correspondence is more subtle. In brief the s-th order large N coefficient 
W±p of the invariant 2r-point functions on a finite lattice can in the compact (+) and 
in the noncompact (— ) model be expressed in terms of a single functional X$?\D](\) of 
the leading order invariant two-point function D - the same functional for both systems. 
This "large N correspondence" is computationally useful because the computations in 
the compact model, which do not require gauge fixing, are much simpler. The 'detour' 
over finite volume cannot be avoided as the correspondence is difficult to interpret di- 
rectly in the infinite volume limit, see [9]. In addition is important to appreciate that 
although the gap equations are related by flipping the sign of the large N coupling 
A = (N+l)//3, the leading order propagator D- in the noncompact model behaves very 
differently as a function of the lattice distance than D + : while the latter shows expo- 
nential decay in the thermodynamic limit, the former decreases only with a power law 
(for d > 2) or increases logarithmically with the distance (for d = 2). 

In a finite volume the gap equations for the compact as well as the noncompact model 
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have the same 0(V) number of solutions for the (mass) 2 parameter; but whereas in 
the compact model there is exactly one positive solution, which is the 'physical' one, 
in addition to multiple negative ones, the sign-flipped gap equation has only negative 
solutions of which only one is the right one for defining D-. The fundamental criterion 
by which the relevant solution is selected should be the stability of the corresponding 
saddle point, but we find also an alternative characterization selecting the physically 
relevant saddle point of the noncompact model. In the thermodynamic limit, in both 
cases all the negative solutions of the gap equation disappear, therefore the large N 
expanded correlation functions, viewed as functions of the lattice points, are not related 
in any simple way, in particular not merely by flipping the sign of A. 

The derivations of the above results also highlight why one cannot expect useful corre- 
spondences to exist beyond asymptotic expansions. For example flipping the sign of the 
coupling in the (exact) generating functionals for invariant correlation functions maps 
the ferromagnet onto the compact antiferro magnet. In contrast, in the perturbative 
asymptotic expansion the sign flip rather relates the ferromagnet to the noncompact 
model, while in the large N expansion the functional dependence on the independently 
defined free propagators gets related. In the latter case the relation to the compact 
model cannot be formulated directly in the thermodynamic limit but one can take term 
by term the thermodynamic limit on both sides of the correspondence. This somewhat 
tricky limit is studied in an accompanying paper [9] in the two-dimensional systems for 
a number of physically interesting invariant quantities. Ultimately we expect the funda- 
mental differences between the compact and the noncompact models to be rooted in two 
facts: the presence of a infinite volume mass gap and the absence of long range order 
(for d < 2) in the compact models and the opposite characteristics in the noncompact 
models for all d > 1. 

The rest of the paper is organized as follows: In section 2 we introduce generating 
functionals for the invariant correlators in the three systems considered and establish the 
perturbative correspondence. In section 3 dual formulations of the generating functions 
are introduced, suited for the large N expansions. The large N correspondence is shown 
first from the generating functionals, then from the Schwinger-Dyson equations, and 
finally verified at low orders. 



2. Invariant correlators and PT correspondence 

Invariant correlators in the noncompact model have to be defined in terms of a gauge 
fixed generating functional. A gauge fixing where one spin is kept fixed turns out to be 
advantageous. To discuss the relation to the compact model we adopt the same gauge 
fixing there. This spin model formulation is also convenient to discuss the relation 
between the exact generating functionals and for the proof of the perturbative (PT) 
correspondence. 



2.1 Definitions 

We begin by setting up the notation and the definitions for the invariant correlation 
functions considered and their generating functionals. We consider the 0(N+1) spherical 
and the SO(l,iV) hyperbolic sigma-models with standard lattice action, defined on a 
hypercubic lattice A C Z d of volume V — |A| = L d . The dynamical variables ("spins") 
will be denoted by n", x G A, a = 0, . . . , N, in both cases, and periodic boundary 
conditions are assumed throughout n x+ L(i = n x . The constraint is n-n = 1 in both cases, 
but with different 'dot' products; namely a-b := a°6° + a 1 6 1 + . . . + a N b N =: a c S c db d in the 
compact model, and a-b := a°b° — a 1 b 1 — ... — a N b N =: a c r/ cd b d in the noncompact model. 
We shall also use the notation a = (a 1 , . . . , o^) for vectors in M. N , so that the bilinear 
forms read a ■ b = a°b° ± a ■ b, in the two cases. Clearly = {n G M. N+1 \ n ■ n — 1} 
is the iV-sphere and H w = {n G M}' N \ n ■ n — 1, n° > 0} is the upper half of the two- 
sheeted N- dimensional hyperboloid. The invariance groups are 0(N + 1) and SOo(l, N), 
respectively. 

Let us briefly note how and are related by symmetric space duality (see for 
instance [10, 11]). Recall that a symmetric space G/K has an involution r associated 
with it such that the Lie algebra q of G decomposes according to Q = t © m as a 
direct sum of vector spaces, where 6 and m are even and odd under r, respectively. 
Furthermore [6, fi] C t, [t, m] C m and [m, m] C t. The dual Lie algebra g* is then 
defined as g* := I © im and the corresponding dual group G* is the (simply connected) 
group whose Lie algebra is q*. G* contains a connected subgroup K having £ as its Lie 
algebra; the dual symmetric space (G/K)* can then be defined as G*/K. On g we have 
an invariant bilinear form B (in our case simply minus the Killing form) which induces 
a dual bilinear form B* on g* . The bilinear forms B and —B*, restricted to m and im, 
respectively, are positive definite in our case and define the metric of the tangent spaces 
at the origins of the symmetric spaces G/K and (G/K)*, respectively. By requiring 
invariance under G or G*, respectively, Riemannian metrics on both symmetric spaces 



are induced. For G = SO(N + 1), K = SO(N),tias gives = SO (l, N) (the universal 
covering group of SO (l, N)) and (G/K)* = SO (l,N)/SO (N) = SO (l, N)/SO (N) 
or G/K = S N and (G/K)* = M N . 



where the upper sign refers to the compact model and the lower sign to the noncompact 
model. In the compact model, with the conventions adopted in (2.3), j3 > models fer- 
romagnetic behavior while (3 < models antiferromagnetic behavior. In the noncompact 
model only (3 > is allowed and the action is unbounded from above, < S-[n] < oo. 



The lattice actions for the two systems are 




(2.1) 



The Laplacian is A xy = — X^[2#x,y — $x,y+p. — 8x,y-fr], as usual. We write 

dtt + (n) = dn N+1 8(n ■ n - 1) , 

dtt^n) = 2dn N+1 5{n ■ n - l)0(n°) , (2.2) 

for the invariant measure on and M. N , respectively. Further 5±(n,n') is the invariant 
point measure on S N , M N , and n T = (1,0, ...,0). Note that the measure d£l + (n) is 
normalized while has infinite volume. 

In the compact model we consider two generating functionals, the usual one and a variant 
with one spin frozen: 

exp W[H] =M j\{ dQ+(n x ) exp { - S + + l - ^ H xy (n x ■ n y - 1)} , (2.3a) 

x x,y 

exp W+ [H] = Af + I Y[ dn + (n x )5 + (n xo ,n1) exp j - S + + i ^ ^(n a • n„ - 1)} {2. 3b) 

where H xy > is a source field and the normalizations A/", A/" + are such that W[0] = 
0. More generally will denote the generating functional with the spin at site 

Xo fixed by a S(n Xo ,n) insertion. In this notation one has = W^t[-H] and 

fdn+(n)expW n [H] = exp W [if]. In fact 



W n [H] = W[H] , (2.4) 

for all n G E> N . This can seen by performing a global rotation n x \— > gn^ =: n x , rr G A, 
with <? G 0(iV + l) chosen such that gn XQ = n\ say. The Boltzmann factor will then 
depend on n x , x ^ x , only and the dQ(n X0 ) = dVt{n X0 ) integration can be performed to 
give I. 

In the noncompact model only the fixed spin variant of the generating functional is well 
defined and we write 

expW^H] = Af_ /jJdfi_(n x )(J_(n a . ,n T ) exp S_ + ]^^H xy (n x - ra y — 1)} ,(2.5) 

x x,y 

where now H xy < sources give damping exponentials. 
Partially connected 2r point functions are defined by 

w±[h] = -j^? w ±A x ^ 2/i ; • • • ; x n Vr) H xiyi . . . u XrVr , 

r>l 

W± :r (x 1 ,y 1 ;...;x r ,y r ):=h xiyi ...h Xryr W±[H] , h xy := — — . (2.6) 

h=o on xy 
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In particular W ±jl (x, y) = (n x ■ n y ) - 1, W ±j2 (x 1 ,y 1 ; x 2 , y 2 ) ■= (n xi ■ n yi n X2 ■ n y2 ) - (n xi ■ 
n yi)( n x2 ' n y2)-> where ( ) are the functional averages with respect to Af± 1 e~ s± . Note 
that W± tr (. . . ; x, x; . . .) = 0. 

In the above we tacitly assumed that W±[.ff] and the correlation functions computed 
from it do not depend on the site xq of the frozen spin and are translation invariant. We 
show now that this indeed the case. If we momentarily indicate the dependence on the 
site as W Xo (and drop the ± subscripts) one has trivially 



Thus, if W xo is independent of x it is also translation invariant. The Boltzmann fac- 
tors (2.3b) and (2.5) can be viewed as a function on the group via F(g , . . . , g s ) = 
f(gon\ . . . , g s n^), where we picked some ordering of the sites Xi, i — 0, 1, . . . , s := V— 1, 
identified n Xi with gin\ and wrote momentarily / for the Boltzmann factor. Then W Xi 
is of the form 



Using the invariance of F under ^ i— > h 1 g i and the unimodularity and invariance of the 
measure dg one verifies that Ui = Uq for all i. 

A peculiarity of the fixed spin gauge is that an invariant 2r-spin correlator (n xi ■ n yi n X2 ■ 
Tly 2 . . . Tl Xr ■ Tly r ), can be re-interpreted as a non-invariant 2r-l-spin correlator. This 
is manifest when Xi = x for some i; since n xo = n 1 * the correlator then involves only 
2r — 1 fluctuating spins. By the above argument any one site can play the role of Xq, 
so that picking such a re-interpretation amounts to picking a site where n Xi = n 1 * . This 
feature holds both in the compact and in the noncompact models, but initially only 
on a finite lattice. When the site x is kept fixed in the interior of the lattice, these 
specific noninvariant 2r — 1 point functions will have a pointwise thermodynamic limit 
in both cases. However the fixed spin averages then do not approach limits which can 
be interpreted in terms of averages without gauge fixing or with a translation invariant 
gauge fixing. This changes if x is identified with a point on the boundary and moves out 
to infinity as A — > Z d . In this case, in dimensions d < 2 an important difference between 
the compact and noncompact models emerges: while in the compact model the Mermin- 
Wagner theorem assures that the non-invariant correlator is equal to its invariant average 
over 0(A^+ 1), in the noncompact model the non-amenability of SO(l, N) prevents this, 
and we find either spontaneous symmetry breaking or divergence of the correlators, as 
discussed in [1]. 

The goal in the following will be to relate the perturbative and the large N expansions 
in the compact and the noncompact models. For the perturbative expansions of the 



W xo [r a H] = W xo+a [H] , (r a H) xy = K 



x+a,y+a ■ 



(2.7) 





(2.8) 



correlation functions, we write 



W ± ^^Y,P~ n wtl. (2.9) 

n>0 

In a large N expansion A := (N+l)/@ is kept fixed and the coefficient functions W^ r in 

^-WTW^^WTir^' (2M) 

are sought. 



2.2 Compact antiferromagnet and noncompact model 

In the literature one finds statements, based on formal manipulations of functional in- 
tegrals, suggesting that compact and noncompact models are essentially related just by 
flipping the sign of (3 (in [2], however, the qualifying remark is added that such a relation 
is 'less clear' at the non-perturbative level). Since this simple 'flip rule' is not completely 
correct and the relations between the compact model, its analytic continuation to (3 < 
and the noncompact model are somewhat subtle, let us try to clarify the situation. To 
avoid confusion, in the following (3 will always be assumed to be non-negative. 

The compact model defined in Eq. (2.1) describes for (3 > a ferromagnet, for (3 < an 
antiferromagnet, so flipping the sign of (3 turns a ferromagnet into an antiferromagnet, 
not the noncompact model. It turns out that one can nevertheless relate the coefficients 
of the perturbation expansion and the large N expansion of the compact ferromagnetic 
and the noncompact model by a sign flip prescription, which does, however, not imply 
any simple relation between the antiferromagnet and the noncompact model. 

On the simple hypercubic lattices we are considering here, there is a transformation 
between the cases (3 > and (3 < 0, based on the fact that we can decompose A into an 
even and an odd sub-lattice 

A = A+UA_, (2.11) 
where x G A + if and only if Yli=i x i ls even > otherwise x G A_. Defining 

e x :=(-l) Xl+ - +a: -, (2.12) 

we can define a map of the configurations by 

n x >-»■ (n' x ) = e x n x , (2.13) 



mapping each configuration of the ferromagnet into one of the antiferromagnet with the 
same (correctly normalized) Boltzmann weight. This yields the following identity relat- 
ing the generating functionals of the antiferromagnetic and the ferromagnetic systems: 



W + ^[H] = W +tfi [H*\ + \ Y^ily - H xy ) + const , (2.14) 



xy 



where we momentarily indicated the dependence on j3 and set H xy := e x H xy e y . For the 
invariant two-point function this gives 



(jl x ' Hy) ^x^yi^x ' fly) 



(2.15) 



and similarly for the higher correlation functions. This implies that the physics of the 
antiferromagnet is essentially the same as for the ferromagnet, in particular they have 
the same mass gap, phase structure etc.. 

The relation between the ferromagnetic 0(iV+l) and the noncompact SO(l, N) systems 
(the latter only exists in the ferromagnetic version) is, however, more subtle: starting 
from the 0(N + 1) model we re-parameterize the variables n x G S w by introducing 
j3 =: g~ x and n x =: ^tt x , so that 



° x n° y + gn x ■ n y = a x a y y/l - gn* ^l-gir$ + gir x ■ n y , 



(2.16) 



where the a x are variables taking the values ±1. To obtain W + [iif] in this parameteri- 
zation one has to (i) integrate over all jf x subject to the constraint 7f^ < /3 and (ii) sum 
over all a x = ±1 except a Xo , which is 1. Explicitly 

exp W + [H] =M + Y, 6^ f J] dn x S(n X0 )9(\P\ - n 2 x ) (1 " 9<)~ 1/2 
x exp < ~y (vx^l-g^ - a x+(lS Jl- gTx 2 x+ ^j + g(n x - jr x+(i ) 
\\ Hx v (^yV 1 -9^1 \j l -9^ 2 y + gn x -n y -lj \ . 

(. x,y ) 



x exp 



(2.17) 



For (3 — > +oo (ferromagnet) the term in which all a x are equal dominates (for (3 — > — oo 
(antiferromagnet) instead the terms with alternating o x dominate). 

On the other hand, we may re-parameterize the noncompact model in a similar way by 
introducing n x = yjg tt x , but this time 



n x ■ n y = n° x n° y - gn x ■ n y = ^l+gn^ Jl+gn 2 - grr x ■ n y . 



(2.18) 



n 



Inserting this into the expression for VF_[if] given in (2.3b), we now obtain 



expW.[H] = N- fl[d7r x 5(7r xo ) (1 + g^)' 1 ' 2 

J X 

X 6XP { Yg ^ (v 71 +^ - V 1 + - - 

te^(^+^v /i+ ^"^'^ _i ) ■ (2 - i9) 



x exp 



Comparing now Eqs (2.17) and (2.19), we see that the partition function of the noncom- 
pact model is obtained from the one of the compact one by (i) flipping the sign of /3, 
(ii) dropping all terms except the one with all a x = 1 and (in) omitting the 9 functions 
restricting the domain of integration. 

So it is clear that the analytic continuation of the compact system to negative j3 (the 
antiferromagnet) is not equivalent to the noncompact SO(l,iV) system, as one might 
infer from a cursory reading of the literature (for instance [2]). This non-equivalence 
will be shown even more manifestly below. Within the framework of the Schwinger- 
Dyson Equations used in section 3.5 a noteworthy consequence is that even the exact 
SD equations do not determine their solution uniquely. 

In the following sections we will see, however, that nevertheless the perturbation expan- 
sions in a finite volume of the ferromagnetic compact and the noncompact model are 
related simply by flipping the sign of (3 (but without the introduction of the sign factor 
e x as in the antiferromagnet); furthermore we will derive a relation between the 1/N 
expansions of the two models. 



2.3 Perturbative correspondence 

We now consider the perturbative expansions (2.9) of the correlation functions. On a 
finite lattice it turns out the perturbative series in the compact ferromagnetic and in 
the noncompact models are asymptotic expansions which are simply related by a sign 
flip in the coupling, (5 \— > —f3; the perturbation expansion of the noncompact model is, 
however, not equal to that of the antiferromagnet, as one might guess from this. 

Concretely 

W W = (-) r+n u>JJ , Vr > 1, n > , (2.20) 

for the coefficients (2.9). The sign flip rule was stated explicitly by Hikami [6] (in a formal 
continuum expansion based on dimensional regularization) but was presumably known 
to other authors in special cases and at low orders, e.g. [5]. It has also been observed 
in the perturbative beta functional of Riemannian sigma models (using dimensional 
regularization and minimal subtraction) and can readily be verified to all loop orders in 
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this framework. But to the best of our knowledge no general proof is available in print, 
so we present a simple proof here. 

We set Z±[H] := expW±[H] and define 



n 



x 



X,fl 



X 



exp {]> H *y (v 7 1 - 9^1 \j l -9^l + g*x ■ T? y - i) \ , 



(2.21) 



and 



Z VT [HW) := /[]^W-a(l + ^)- 1/2 

J X 

V x,n 

\ \ Yl H xy(V 1 + 9 j ?l \j l + 9K 2 y - g^x ■ T? y - lH . 



x exp 



(2.22) 



Z^ T [iJ] is, up to a multiplicative constant, just the term of Z + [H] (Eq.(2.17)) in which 
a x = 1, for all x. Z^ T [i/] differs from Z_[H] in Eq. (2.19) by the presence of the 6 
functions and a multiplicative constant. We now state the 



Result: 



Z ± [H]=c ± ((3)Zr[H}(P) + 0(P 



(2.23) 



where c±(f3) are normalization constants independent of H . Further the Z^[H](f3) have 
asymptotic expansions of the form 



Zi T [H]((3)~J2an[H](±(3)-\ 



(2.24) 



n=0 



with identical coefficients a n [H] in both cases. 

The existence of the asymptotic expansions (2.24) is a standard result of Laplace's 
method in asymptotic analysis (see for instance [12, 13]). The relation (2.23) expresses 
the irrelevance of boundaries away from the maximum as well as of the terms in which 
the o x are not all equal. It follows from two facts: 

Fact 1: On a finite lattice in the compact ferromagnetic model all terms in Z + [H] except 
the one with a x = 1 are 0(f3~°°). 
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Fact 2: On a finite lattice in the noncompact model for any x G A and (3 > (3 > the 
expectation value 

(9{n 2 z -l)) := j\\d^{n x )5{n xo ^)0{nl-l) 

L J X 

x exp | -S- + ]-^H xy (n x - n^-l)}, (2.25) 
is bounded by aexp(— bf3) with a, 6 > . 

Corollary: There are constants a±, b± > such that for (3 > (3 



ktV-1 



Z™[H](P)p- N . 



Z±[H]((3) 



< a± exp(-6±/3) . (2.26) 



Both facts have essentially the same origin: the fact that the 'energy' in the omitted 
contributions is exponentially small. We give a detailed proof only of Fact 2 and add 
some comments on the proof of Fact 1. 

Proof of Fact 2: We do the integrations in the following way: first we integrate over all 
configurations having a fixed value s of the action and afterwards over s. Thus we can 
write 

Z-[H\{(3) = Jdse- ps PH {s) , (2.27) 

with some non- negative density function p#. If we then can show that for any config- 
uration violating the field cutoff the action is larger than s ('energy bound'), it will 
follow that 

Z.[H}((3) (6(n 2 x - 1)) < J dse- /3s p H (s)e(s - s ) . (2.28) 

The right hand side is easily seen to be < exp(-fso) x Z_[H\((3/2). By standard 
asymptotic analysis (see for instance [12, 13]) it follows that Z_[H](f3) /Z_[H](f3 /2) goes 
to a finite limit as f3 — * oo, and therefore 

Z-[H\{(3) - a ' {2J9} 
for f3 > fy. So Eq. (2.28) implies Fact 2. 

It remains to prove the 'energy bound'. Since we are using the fixed spin gauge, we 
have n xo = n) . Then by assumption there is a site x\ with n Xi > \/2. Choose a path 
P of length |P| < V from xq to x\ and denote the part of the action corresponding to 
the path P by sp. We claim that sp is minimized by moving in equal steps along the 
geodesic from n XQ = to n Xl , which implies 

8 > s P > > ^ , (2.30) 



1 1 



with some constant c. 

To see the lower bound on sp it is sufficient to consider a path of three points, since the 
minimization condition is local. Thus we only have to minimize e(n) := riQ-n + n-ni over 
n and show that the minimum is assumed for n being the midpoint of the geodesic in 
M N from no to n\. Choosing coordinates such that n = w) and n\ = (ch#i, sh#i, 0, . . . 0) 
with 6*i > we find 

e(n) = n°(l + ch^) - n 1 sh6 l > n°(l + ch^) - ^(n ) 2 - 1 sh0i . (2.31) 
Putting n° = ch.9 this becomes 

e(n) > ch9 + ch(9 1 - 6) > 2ch^- , (2.32) 
using the convexity of ch function in the last step. It follows that 

S p>|P|(ch(V|P|)-l)>|L>^, (2.33) 

with c = (arccoshv^2) 2 /2. The proof of Fact 2 is now complete with a as above and 
b = c/V. 

To show Fact 1 one proceeds essentially in the same way; only the energy bound is more 
easily derived in a slightly different way: again it suffices to consider three spins no, n, n\ 
with the energy 

e(n) = n ■ no + n ■ n\ , (2.34) 
direct minimization over n yields 

n = A(n + ni), (2.35) 

with A fixed by the requirement n 2 = 1. Son has to lie on the great circle connecting n 
and ri\. As before we conclude from this that the minimal energy for a path is obtained 
by moving along the shorter part of a great circle (geodesic) connecting n to rt\ in equal 
steps. 

The corollary as well as the main result are obvious consequences. 

Thus we have learned that the perturbative expansions of the compact ferromagnetic 
model and the noncompact model are related just by flipping the sign of f3, even though 
doing the same sign flip to the full model leads to a completely different model - the 
ant if erromagnet . 



1 o 



The relation between the asymptotic expansions in 1/(5 for the compact and noncompact 
models is given in (2.20). For the two point function this means for instance 

This should be contrasted with the relations between the perturbation expansions of the 
two-point function between the ferromagnetic and anti ferromagnetic compact models, 
which follows from (2.13) 

W-\{X, y) lantifcrro = e x 6 y W^\(x, y) |f C rro • (2.37) 

It is instructive to see what this means for the lowest order (tree graph) asymptotics of 
the two point function: 

(n x -riy) = 1 + -w^ tl (x, y) + 0(/3' 2 ) ferromagnet 

(n x ■ n y ) = ^1 + -^w+^x, y)^j e x e y + 0(j3~ 2 ) antiferromagnet (2.38) 

(n x ■ n y ) = 1 — — w+^x, y) + 0{(3~ 2 ) noncompact model , 



where 



with E p as in Eq. (3.22) below. 

Note that although VF + [if] for the ferromagnet and the antiferromagnet are related by 
flipping the sign of j3, this is no longer true for the asymptotic expansions, even in 
finite volume. Rather the sign flipped perturbative expansion is that of the noncompact 
model. 



3. Dual formulations and large N correspondence 

The goal of this section is to establish a relation between the large iV expansions of 
invariant correlators in the compact and the noncompact sigma-models. In the compact 
model the large N expansion is based on the familiar dual formulation, which here has to 
be modified due to the fixed spin. In the noncompact model the duality transformation 
is ill defined but by performing Gaussian integrations in horospherical coordinates, one 
can obtain a well-defined counterpart of the dual formulation. 
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3.1 Noncompact generating functional via horospherical coordinates 



Before turning to the expansions proper we present here an exact rewriting of 
based on a partial evaluation where the (dimHAr) y integrations are reduced to V integra- 
tions. The point of departure is the fact that the hyperboloid M N admits an alternative 
parameterization in terms of so-called horospherical coordinates. These arise naturally 
from the Iwasawa decomposition of SOo(l,iV). Here it suffices to note the relation to 
the hyperbolic spins 



N-1 „ N-1 



n ° = c h6 + - *i e ~* > n 1 = ak6 + -^2^e- e , n* = e~ Vi, i = 2, . . . , N , (3.1) 

i=i i=i 

and that (9, t±, . . . , ijv-i) & defines a globally valid system of coordinates. It is 
convenient to write t — (t±, . . . , ijv-i) and t-t 1 = tit' x + . . . ijv-i*jv-i- F° r the dot product 
of two spins n x , n y G M N this gives 

n x ■ n y = ch(e x - 9 y ) + \{t x ~ t y ) 2 e- e *- e y , (3.2) 

and for the measure in (2.2) 

dn_(n) = e-W-VdOdt! . ..dt N -! = e-W-VdOdt. (3.3) 

The key advantage of horospherical coordinates is manifest from (3.1), (3.2): for a 
quadratic action of the form S- in (2.1) the integrations over the t variables are Gaussian 
and can be performed without approximations. We refer to [4] for the derivation and 
only note the result 



expW_[#]=exp{- ^Vtf^W I T\da x 

1 2 x,y } MH) x ^ xq 

x exp { - ^Trln A - t £ a x + ^{A-Xl } ■ (3.4) 

X^XQ 



Here 



^■xy ^xy rflxy &xy&x A X y -\- Cl X0 S XX0 S X y , (3-5) 



P 



and A xy is the matrix obtained from A xy by omitting the xo-th row and column. The 
domain T>(H) is an algebraic variety described by 

V{H) = {ae [-2d, oo] v ~ l | det A > , A positive semidefinite} . (3.6) 
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We also anticipate from [4] the following 

Result: The correlation functions W_^ r admit an asymptotic expansion of the form 
(2.10), whose (uniquely defined) expansion coefficients coincide with those defined by 
the Laplace expansion of (3.4) where T>(H) has been replaced by R^ -1 . In turn these 
coefficients coincide with those of the formal large N expansion of 

expWl[H] = exp { - ^J2 H w} ^ f T[ da * ex P { _ (N+l)S-[a, H] } , 

x,y x^xq 

X^XQ 

where a x corresponds to 2i\a x . However 

W-[H]^W^[H]. (3.8) 

Heuristically the origin of the 'dual' generating functional W^if] can be understood by a 
dualization of the 'spatial' spin components n x , x G A, and a formal contour deformation 
in field space [9]. The formal nature of the latter is responsible for (3.8), although the 
large iV expansion coefficients of W- [H] are correctly reproduced. 



3.2 Dual formulation: compact model 

In the compact model the counterpart of (3.7) is obtained along the familiar lines: 
one first implements the constraints n x -n x = 1 via a Lagrange multiplier field and 
then performs the Gaussian integrations. Of course in the compact model no gauge- 
fixing is required and the result is well-known. For the purposes of comparing with 
the noncompact model, however, we want to perform the dualization here in the fixed 
spin gauge, i.e. for in which case certain modifications occur. The resulting dual 

generating functional W^[if] is an exact rewriting of the original one: W[if] = VF+[if] = 

iv;'[//]. 
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In the fixed spin gauge Gaussian integrals of the following form arise 
/ II d(j> x 6((f> X0 ) exp | - - ^2 <f>xA xy (j>y + ^2 J x<t>x\ 

X x,y X 

= (27r) Z ^(detl)- 1 / 2 cxp{i^J x (A- 1 )^J J/ }, (3.9a) 

x,y 

d(j) x exp | - - ^ 4>xA xy 4> y ^ 

x^=x x,y 

= (2 7 r)V(detl)-^exp{- ^' }, (3.9b) 

for a real field <p x , x G A, and a symmetric invertible L d x L d matrix A, for which we 
assume that (A~ V ) XQXQ ^ 0. Here A is the (L d — 1) x (L d — 1) matrix arising from A by 
deleting its xo-th row and column. The inverse of A can be expressed in terms of the 
inverse of A via 

(A-\ y = (A~% y - i^^lm (3 10) 

ysi )x x 

This equation makes sense a priori for x, y, ^ Xq, but we may trivially extend the matrix 
A^ 1 to a L d x L d matrix by setting ([A^ 1 ]^)^ = in accordance with Eq. (3.10). Then 
A and [yl _1 ] cxt are of course not inverse to each other 

^ ] A xz ( [A ] C xt)zy = $xy ~ $xxo ^-1^ ' (^-H) 



XqXq 



To verify (3.9) one writes Ex.j,^^^!/ = Ei^M^j + ^loEj^^o^i + 
^xo^^o^o- To get (3.9a) one first does the trivial <j) Xo integration. To get (3.9b) one 
uses (3.11). The determinant of A is related to that of A by 

det A 

tetA = -— . (3.12) 

ysi jxoxo 

Often a term in the x -th matrix element on the diagonal of A has to be split off according 
to A xy = A xy — c5 xy 5 XQX . In this case the inverse of A is related to the inverse of A by 

(A l ) xy = (A l ) xy + - ~ ^xxoiA 1 ) ?/a ;o • (3.13) 



In particular A xoxo - (A x ) 1 = A xoxo - (A *) 1 and 



. ^} xxo = ( ~ 1)3:3:0 . (3.14) 



XqXq 



With these preparations at hand the dualization of VF + [if] is straightforward. It is 
instructive to start from single out one spin, n Xo , and to postpone its dualization. 

Inserting 

/7 
exp{-t{N + \)a x [n\ - 1)} for x ± x , (3.15) 

the interchange in the order of integrations can be justified and the relevant Gaussian is 
of the form (3.9b) 

/n*.«p{-§^£^} 

x^£x x,y 

= Const (det A)-^ exp { - ^^+1^)^} , 

A xy = -A xy + 2i\a x (l - S XQX )S xy - j^—^H xy . (3.16) 
This gives for the dual generating functionals 

e W W d [H] = Jdn Xo 5(nl o -l)e W W* o [H], 

exp< o [tf] = exp{-i^^}Ar|n^ ex p{-( Ar + 1 )^„[«^]}' 

x,y x^xq 

S^ [a,H] = iTrlnA-^^ + ^^I- 1 )- 1 ^. (3.17) 

X^XQ 

Here we used the fact A = A; [H] is the generating functional in the fixed spin 
gauge and [cx,H] is its dual action. We shall also write S + [a,H] for S n i[a,H] and 
W^fiJ] for W^ T [i?]. In contrast to the noncompact model one has 

W*[H] = W + [H] , W d [H] = W[H] , (3.18) 

and also W+[H] =W[H], by (2.4). 
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As a check on (3.17) one can verify that by dualizing also the last spin one recovers the 
familiar expressions. Indeed, 



dn xo 5(n 



XQ 



-)exp{-^(I-\-U<} 



/ 2ttA \(n-i)/2 r 

= \W+i) J 



2A 

da Xn e i{N+1)a ^ 



(A- 1 )-'+2i\a 



X'O 



-(AT+l)/2 



(3.19) 



Using now (3.13) for c = —2i\a xo and (3.12) one obtains 

N + l, 



expW[H] = J\f jY[da x exp{ - N * l r Tr In A + i ^ a a | , 



= -A x „ + 2^0^ 



A 

ivTT 



(3.20) 



as required. 



3.3 Finite volume mass gap and basic propagators 

The saddle point expansions of W± [H] define to leading order the gap equations and the 
invariant two-point functions. In both quantities the dependence on the site xo of the 
fixed spin drops out, see section 3.6. The result for the leading order two-point function 
can be written in the form 

XWi%,y) = ±\D ± (x -y)-l = {n x - n y )\ N=oo - 1 , (3.21) 

where D±(x — y) is the basic propagator. For the compact model it has the well-known 
structure D + {x) = D(x)\ ul=UJ+ , where 

1 «^ p ip ' x 

with the sum over all p = ^-(^i, • • • , no), rii — 0,1, ... , L—l, and E p := 2d — 2 cos(p • 
fi). Further uo + = uj+(\, V) is the dynamically generated finite volume mass 2 term 
determined by the gap equation -D(O) = 1/A. A subtlety is due to the fact that the 
equation for cj+(A, V) is an algebraic equation of degree 0(V) and therefore has 0(V) 
solutions. Among those solutions there is always a unique positive one, but it is not 
obvious that this is always the physically relevant one: in dimension d > 3 for weak 
coupling there is spontaneous symmetry breaking and it has been found in [14, 15] that 
at least in the determination of the constrained effective potential one has to take a 
negative solution of the gap equation. Here we want to eliminate these complications 
by always working in the 'strong coupling' regime, which we will now specify. 
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In the infinite volume limit the gap equation becomes 



and it has no negative solutions. It has a unique positive solution for A > \*(d), where 
A* (2) = and X*(d) > for d > 3 (A*(d) is known as the critical coupling of the spherical 
model). From now on we will assume A > A*(d) and always require u + (X,V) > 0. In 
this case it is possible to show that the finite volume corrections u + (X, V) — u + (X, oo) 
are exponentially small. 

In the noncompact model one has similarly D_(x) = D(x) \ L0=UJ _ , where u;_ is a solution of 
the gap equation XD(0) = — 1, see Eqs. (3.54), (3.56) in section 3.6. Clearly analytically 
continuing the multivalued algebraic function uj + (— X,V) to negative values of A and 
setting U-(X, V) := X, V) will produce solutions of this gap equation. But again 
there are 0(V) such solutions, all of which are negative. To see this it is convenient to 
split off the zero momentum mode and to write 

0^ne[0,L-l]d 

where u n = E2-nn.lL , n = (ni, ■ ■ ■ , n<i)- Here f(ou) is a meromorphic functions with simple 
poles along the negative real axis. The first pole is at u> = — 4 sin 2 (7r/L), the others are 
at approximately integer multiples thereof. Specifically 

u n = -A{ii/L) 2 {n\ + ...+n 2 d )+ 0(1/ L 4 ) . (3.25) 

In terms of the function / the gap equations ±XD(0) = 1 read 

Xf(w) = ±1 - -A- . (3.26) 

In the compact model (upper sign) is has a unique solution u + (X,L) > for any given 
A > A*(d); all the other 0(V) solutions are negative. In the noncompact case there are 
only 0(V) negative solutions, one in each interval of the ui axis between consecutive 
poles of /. 

To obtain an asymptotic expansion in 1/N one has to pick the solution corresponding 
to the absolute maximum of the dual Boltzmann factor. But there is a simpler way 
to do this selection: if the large N expansion is to yield an asymptotic expansion of 
the invariant two-point function, we have to demand that — XD(x) > 1 for all x. Since 
(n x -n y ) > 1 in the noncompact model, this is a necessary condition for the interpretation 
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of D_(x — y) = D(x — y)\ u =w_ as the leading order two-point function. This condition 
fixes the solution to lie in the interval u G (—4 sin 2 n/L, 0). Indeed, — XD(x) > 1 implies 

^2(l-cosq-x)(-XD(x)) > V, q^O, 

X 

11 1 , 

> T , q^O. (3.27) 



-uVl+u/E q ~ X 

Since uo < both factors on the left hand side must be positive. Hence E q + uo > for 
all q 7^ 0, for which uj > —4 sin 2 n / L is a necessary and sufficient condition. 

It is easy to sharpen this bound to 

w > "mTT™ 2 ! (3 - 28) 

by noting that 



> 



> "F > - + FT • 3 - 29 



A E p + u uj 4 sin j + to 



From now on we write tu_(X, V) G (—4 sin 2 n/L, 0) for the unique root of the XD(0) = — 1 
gap equation such that — XD_(x) > 1 for all x. Its large volume asymptotics comes out 
as 

i^(l + 0(l/lnV)) d = 2, 

Vu4\V) = { (3.30) 



Q + Cd) \0(V-^) d>3. 



where C7, = J* 

In summary, to leading order the invariant spin two-point functions in the compact 
ferromagnetic and in the noncompact model can be written in the form (3.21), where 

D±(x) := £>(x)|u,=u,±(A,v) , with uj±(X, V) solution of XD(0) = ±1 . (3.31) 

In both the compact model and the noncompact model there is a unique root of the 
gap equation XD(0) = ±1 such that XD±(x) < ±1 for all x. From now on the symbols 
D±(x) will always refer to the unique propagators with this property. They satisfy 

1 < -XD_(x) < -1 - 1X 



2A 

-1 + - — <XD + (x)<\. (3.32) 

The uj± dependent bounds follows from D(x) > 2/(Vuj) — D(0). In particular < 
-cj_ < X/V. 



on 



For V — > oo, -D+(^c) approaches the massive free propagator of squared mass u + (X, oo), 
while D_(x) becomes massless and naively appears to be ill-defined. If one inserts, 
however, for uj the L-dependent solution of the gap equation — A.D_(0) = 1 one can 
rewrite \D_(x) in such a way that the existence of the thermodynamic limit is manifest: 



n„ 



= -XD. 



[x 



N=OD 



y) 



A x ^ 1 



cosk-(x — y) 



(3.33) 



which for V — > oo has the well-defined limit 



1 + A 



/ 



d d k 1 — cos A; -(a; — y) 



(3.34) 



3.4 Large N correspondence 

The results of sections 2 and 3 now lead to the 

Result: Let W r>+ and W r - denote the connected parts of the invariant 2r-point func- 
tions (n xi -n yi . . . n Xr -n yr ) in the SO(A^+ 1) and in the SO(l, N) nonlinear sigma-model, 
respectively. (With bare coupling f3, defined on a hypercubic lattice A C Z d of volume 
L d with periodic boundary conditions, and a suitable gauge fixing in the noncompact 
case.) Then the W± jr have well-defined "large N" asymptotic expansions of the form 
(2.10), with A := (n'+ l)/(3 fixed. Further: 

(a) The coefficient functions W±^ r are translation invariant. The lowest order two-point 
functions have the form W±\(x,y) = ±D±(x,y) — l/\, where D±(x) = D(x)\ U! ^ ul± . 
Here D(x) is the free propagator of squared mass uj and oj±(X, V) are the unique 
solutions of the gap equations ±\D(0) = 1 with the property ±\D(x) > 1, for all 
x. 

(b) For all r > 1, s > 0, there exist unique functionals Xr- S \D](X) of a free propagator 
D such that W+\ = Xr S \D + ](X) are the coefficients in the compact model and 
W[ s) r = (-l) r X r (s) [,D_](-A) are the coefficients in the noncompact model. 

Proof. The existence of the large N asymptotic expansion has been shown in the compact 
model in [16], for W[H\ (or its the counterpart with a linear source coupling). Since 
W^[if] = VF d [iT] = W[if], this directly carries over to VF+[if]. In the noncompact 
model the result has been anticipated in Subsection 3.1; for the derivation we refer to the 
forthcoming paper [4]. The two statements in part (a) follow from sections 2.1 and 3.3, 
respectively. For (b) we use the result stated in section 3.1, namely that W d [H] (although 
not equivalent to W_ [H] and not obtainable from a well-defined duality transformation) 
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generates the correct large N asymptotic expansion. Given this, it suffices to note that 
in (3.17) and W^if] in (3.7) have the same structure with actions that are 
related by the involution 



Hence, if indeed the moments of [if] and [H] have asymptotic expansions of the 
form (2.10), their coefficients are uniquely defined and must therefore also be related by 
the involution (3.35). 

For the sake of clarity let us add two remarks: 

First, VF + [if] and do not coincide exactly, as follows from (3.8), (3.18) and the 

fact that the boundary T>(H) in (3.6) is invariant under the involution (3.35). In physics 
terms this amounts to the observation made before that the antiferromagnet and the 
noncompact model are inequivalent. 

Second, (b) does not state that the large N coefficients of the compact (ferromagnetic or 
antiferromagnetic) model and the noncompact model are related by flipping the sign of 
A. The latter is certainly incorrect. Rather the dependence on D± is same after flipping 
the sign of A, but D + and D_ are very different and independently defined functions 
of the lattice distance. So, given the large N coefficients of W +>ri as a function of the 
lattice points, there would be no simple way to construct the large N expansion of W-^ r . 



3.5 Schwinger-Dyson equations 

It is instructive to look at the previous result from the viewpoint of the Schwinger-Dyson 
(SD) equations. Under the assumption that asymptotic expansions of W- r and W +>r 
are known to exist, this also provides an alternative derivation of the result. 

In the compact model the Schwinger-Dyson equations for the moments of VF[if] have, to 
our knowledge, first been formulated by M. Liischer [17]. In the noncompact model the 
necessity to gauge fix requires modifications. In the fixed spin gauge the modifications 
in the equations are minimal, but they inevitably refer to the preferred point xq. In 
the following we formulate the SD equations in parallel for the gauge fixed functionals 



S + [a, H] i— > S-[a, H] for a x h- > — a x , A h- > —A . 



(3.35) 



W ± [H}. 



The result is: 



±(3A Z [h zy W ± - h zx W ± - h zx h xy W ± - (h zx W ± )(h xy W ± )\ | 



\z=x 



+ Yl H ** h *v W ± - h zxW± - h xy W± - h xz h xy W± - (h zx W±)(h xy W±)\ 



-N(l-5 xy )(h xy W± + l) =0. 



(3.36) 



Again the upper sign refers to the compact model and the lower one to the noncom- 
pact model. In addition there is the following minor but essential modification: in the 
compact case (3.36) holds for all 1,1/ 6 A, since W+[-ff] coincides identically with the 
generating functional defined without gauge fixing, see (2.4). In contrast, in the 

noncompact case it holds for all x$ 7^ x G A and all y G A. Specializing to y = xo thus 
results in an equation that is qualitatively different from the others. 

The derivation of (3.36) is based on the invariance of the product measure Y[ x dQ(n x ) 
under rotation of any one of the spins, n xi , say, where X\ 7^ x in the noncompact model. 
For an infinitesimal rotation with the Lie algebra element t ab we write 

6%ni = Sunlit*)*. (3.37) 

For later use we pick an explicit basis and note the completeness relations: 
so(N + 1) : (t ab ) d = 5 a c 5 bd - 5 c b 5 ad = (t ab ) d , 

\ Y,( tab )c d (t ab )J = M /d - W . (3.38a) 

ab 

80(1, N) : (t ab ) c d = 5 a cV bd - 5 c b V ad = V aa 'v bb '(t a , b/ ) c d , 

\ J2( tah )c d (tab) e f = VceV fd - S c f S e d , (3.38b) 

ab 

where 77 is defined by a c i] cd b d = a°b° — a^b 1 — . . . a N b N . Denoting by ( )# the source 
dependent averages this gives in a first step 

(-(8?S)0 + J2 H xz t ab n x -n z O + 5 ab O) H = , (3.39) 



for any local function O = 0({n}) of the spins, and x 7^ xq in the noncompact case. 
Specializing to O — n x ■ t ab n y and summing over a, b, using the completeness relations 
(3.38) gives 

(± (3A z [n z ■ n y - n z ■ n x n y ■ n x ] + V] H xz [n z ■ n y - n z ■ n x n y ■ n x \ 

\ z=x 

- N(l - 5 xy )n y ■ n x ) =0. (3.40) 

Replacing the ( )h averages by functional derivatives of W±[-ff] results in (3.36). 

Coupled bilinear equations for the W r functions arise by differentiating (3.36) repeatedly 
and setting H equal to zero. The equation obtained by s-fold differentiation involves all 
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W r functions with r < s + 2; as a consequence no finite subsystem of equations arises 
for any individual W r . 

Before turning to the large N expansion of the system (3.36) we wish to stress that 
even the exact SD equations do not determine their solution uniquely. To see this, let us 
momentarily denote the functional equations (3.36) with the upper and the lower sign by 
(SD)+ and (SD)_, respectively. Further we write W+,p, W+-p, and W- : p, (3 > 0, for the 
generating functionals of the ferro magnet, the anti-ferromagnet, and the noncompact 
model, respectively. Then, by construction, W +t p solves (SD) + while both W +t -p and 
W-^ solve (SD)_. (Consistency requires that substituting the rhs of (2.14) into the 
(SD)_ equation converts it back into a (SD) + equation for W +t @ evaluated on H e . Using 
(Aef) x = —e x (Af) x — Ade x f x , this can be verified to be the case.) Thus if the exact 
SD equations were to determine their solution uniquely, W+-p and W- t p would have to 
coincide. This however contradicts the discussion in section 2.2, viz the antiferromagnet 
and the noncompact model are physically and mathematically inequivalent. The upshot 
is that somehow initial or boundary conditions have to be imposed on (3.36) to specify a 
solution uniquely. Since the exact equations couple all multipoint functions, this seems 
difficult to do concretely. 

In contrast, the large N ansatz (2.10) effectively converts the equations into ones which 
can be solved recursively, and 'initial' conditions can be specified. The structure of the 
large N expanded SD equations is best explained by spelling out the first few (see [9]) 
from which one can read off the recursion pattern for the Wr S \ r + s > 1, functions: 

Wf> ->W< 0) ->Ws (0) -> 
I I 

i 



Fig. 1: Recursion pattern for the solution of the large TV expanded SD equations. 



To compute a given coefficient all quantities having arrows pointing towards it are 
needed. The detailed form of the equations and the solutions is not essential for the 
following argument. 

We only need: (i) the fact that the equation for is autonomous and reads 



±A z [wi°\z,y) - W?\z,x)\ =F XW^(x,y)A g W^(z,x) 



(0), 



r(0), 



l-5 xy + XWi°\x,y), 



(3.41) 



and (ii) the assumption that each recursion step in Fig. 1 has a unique solution. We 
first present the solutions of (3.41) and then discuss the status of the assumption (ii). 
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A solution of (3.41) with the upper sign is 



Wi°\x, y) = D + (x -y)- 1/A ferromagnet , (3.42) 

Two solutions of (3.41) with the lower sign are 

Wi°\x,y) = e x e y D + (x — y) — 1/A anti-ferromagnet , 

Wi°\x,y) = —D_{x — y) — l/\ noncompact model , (3.43) 

where for the verification of the first solution (Aef) x = —e x (Af) x — Ade x f x has been 
used. For each solution (x, x) = amounts to the corresponding gap equation. As 
discussed in section 3.3 the gap equations have 0(V) solutions other than the physical 
ones u± entering D±. Each solution of the gap equation defines a different solution 
of (3.41) via (3.42), (3.43), and in principle there could be others. The selection of a 
specific solution - here D± - can be justified either from the functional integral (in that 
it gives rise to valid asymptotic expansion) or from the physics one seeks to describe. 

The uniqueness of the recursion (ii), to all orders, is probably difficult to establish 
directly from the SD equations. In a perturbative expansion of the W r the problem 
becomes linear and the fact that the perturbative expansions of the W r is uniquely 
determined by the equations has been pointed out by M. Liischer [17]. In other words a 
potential nonuniqueness is known to be of order 0([3~°°), directly from the SD equations. 

An indirect way to establish the uniqueness (ii) is by showing that the moments of 
the generating functionals VF±[if] have well-defined asymptotic expansions of the form 
(2.10). The expansion coefficients then must be unique, and since W^fif] are solutions 
of (3.36) the coefficients W±\ must be solutions of the expanded SD equations, and 
the unique solutions. This also implies that the solutions must be translation invariant, 
despite the preferred role of the x = Xq equation. As mentioned earlier, the fact that 
the moments of W + [if] have an asymptotic expansion follows from [16]; for W^-fif] this 
will be shown in [4] 

Given the uniqueness of the recursion one can readily re-derive the result of section 3.4. 
From (3.36) and (2.10) it is clear that the involution 

A h-> -A , W^ n \\) h-> (-) r W^ n \-\) , (3.44) 

will map the Schwinger-Dyson equations of the compact model onto those of the non- 
compact model, to all orders of their large N expansions. To draw conclusions about 
the solutions, however, the uniqueness asserted in (ii) is crucial. It implies that a se- 
quence of solutions Wr n \\) in the compact model, based on a specific choice of W^l(X), 
has a unique counterpart in the noncompact model based on the flipped initial solu- 
tion — A) of the lower sign SD equations (3.41). Since the physical solutions 



o rr 



of the compact ferromagnet and the noncompact model are related in this way it fol- 
lows that there exists unique functionals Xr S \D](X) such that W+ \ = Xr S \D + ](\) and 

W ( _ s l = (-l) r xi s) [D_](-A), where D±(x - y) := ±{wf\{x, y) + 1/A) are the physical 
solutions of (3.41). In summary, the previous discussion reproduces the parts (a),(b) of 
the result in section 3.4. 

We add some remarks. First, the involution (3.44) is the counterpart of (3.35). Let us 
repeat that the involutions (3.35) or (3.44) do not imply a simple relation between the 
Wr n \\) and Wr n \— A) as functions of the lattice points, because already Wi°\\) and 
are based on solutions of different gap equations. It is rather the functional 
dependence on D± which is the same. 

The involution also does not map the large N coefficients of the compact model into 
those of the antiferro magnet, defined by Vy + [iJ]|^_>_^. The reason is that the large N 
saddle point relevant for the antiferromagnet is not the one obtained by reversing the 
sign of A in the gap equation. Rather is is given by the same equation as the saddle 
point for the ferromagnet, as is manifest from (3.43). Generally the large N coefficients 
of the antiferromagnet are obtained simply by introducing appropriate products of the 
alternating sign function e x in those of the ferromagnet, as explained in Section 2.2 . 

The special role of the SD equations (3.36) for y = xq is related to the discussion in 
the paragraph following Eq. (2.8). It may suffice to illustrate the point with the lowest 
order equation (3.41). The solution of the y = xo equation in (3.41) is 

("S)/.-. = «x + 0(^ I ) with n x :=l + XWi°\x ,x). (3.45) 

The interpretation of this solution as the nonzero average of the n° x component follows 
from [1]. To higher orders the solutions of the y = xq equations produce corrections to 
(3.45) or non invariant three-point functions, etc. 



3.6 Algorithm in fixed spin gauge 

Here we illustrate the main results from sections 3.1 (existence of the asymptotic expan- 
sion) and 3.4 (compact-non-compact correspondence) by an explicit computation of the 
two- and four-point functions in the noncompact model to next-to-leading order. The 
starting point is the generating functional (3.5), where according to the result described, 
T>(H) has been replaced by Due to the gauge fixing a number of complications 

arise compared to the conventional large N computations; in particular it is not obvious 
how translation invariant correlation functions are recovered. 

Initially we keep the putative saddle point configuration u x ,x ^ x , generic and spe- 



cialize to u x = cj„ only later. We thus write 

a x =u x + | > G K , a: 7^ a; . (3.46) 

Although notationally cumbersome it is important to carefully distinguish "hatted", 
"tilde" and "plain" matrices. We set 

M X y . A-xy ^x$xy • M X y + ^xo^xxo^xy i D X y . )xy ■ 

Uxy ■ u x S X y • U X y ~\~ u Xo 5 XXQ 5 X y . (3.47) 

The matrices M and U are also defined for x, y — x . For M we assume u Xo ^ uj- + A, 
which ensures its invertibility; in U xy the variable u Xo is non-dynamical. Further we 
write M, U for the matrices obtained from M, U, respectively, by deleting the XQ-th row 
and column and put D = M~ x . Then 

A=M+ - 1 U+—^—\H, 
v/iVTl N + 1 

A~ 1 =(l + . 1 du+—}—XDh) 1 D, (3.48a) 

V v^vTT N + i ) v 1 

A=M+ - 1 Q+-J—XH, 
y/WTT N + 1 

In A=\n M + In ( 1 + 1 n/7 + — - — \DH^\ . (3.48b) 
The required expansions are for square matrices a, b 

i 1 + 7FTT + " = 1 + g WTTF * (a ' 6) ' 

5i(a, b) = — a, q 2 (a,b) = a 2 — b , g 3 (a, 6) = — a 3 + a6 + 6a , 

g 4 (a, 6) = a 4 - a 2 6 - a&a - 6a 2 + o 2 , (3.49a) 

1 11 

pi(a,6)=a, p 2 (a,b) = --a 2 + b, p 3 (a,b) = -a 3 - -(ab + ba) , 

p 4 (a, 6) = -^a 4 + -(6a 2 + a&a + a 2 o) - -o 2 . (3.49b) 
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Explicit formulas for the polynomials qi,pi, I > 5, can easily be found but will not be 
needed. With this notation one has 



(a- 1 ; 



xqxq — D XqXo 



/> 



(3.50a) 



Qi:=D; 2 X0 { qi (DU,\DH)D) 



XqXq 



ln 1 = ln M + E {N +iy,2 Vi (DU, XDH) , 
from which the expansion of S[a, H], i.e. 

SW,H] = ^lnA+^^a,-^^- 1 )- 1 ^, 



= {N 1 1)l/2 Si[u, 



m 



(3.50b) 



(3.51) 



can readily be computed to any desired order. We present Si, . . . , S4 first in the con- 
densed form arising from (3.49), (3.50). 



S 1 [u, H] = Itt Pi 0U, XDH) + ±.J2 u * + T^Qi > 



(3.52) 



2A V 2A 

5 2 [«, H] = ^Tr P2 (DU, XDH) + -L(Q 2 - A^Q 2 ) , 
5 3 [«, #] = ^ P3 (5C/, AM) + ^(Q 3 - D X0XO 2Q 2 Q 1 + D 2 X0X0 Ql) , 

5 4 [«, HI = ^Trp 4 (5f?, M) + ^(Q 4 - 5 X0X0 (2Q 3 Qi + Ql) + 35 2 o:co Q 2 Q 2 - 5jU<#). 



Note that 5; is a homogeneous polynomial of order / in the u x s and fTs, assigning tt a 
degree 1 and H xy degree 2. 

For Si and 5 2 this gives in a first step 



2 i 



(3.53) 



S 2 h#] = 

x,y^x 



^x^y 



\f) 2 -I- ' Dxpx j-. Dx y 



1 ~ D 2 D 2 



2\d y D 2\ M)M) d 2 D 2 

^XQXQ ^XQXQ ^XOXO ^ X0X0 



+ H xy 

x,y 



\ n D xox D xoy 



D X q X0 D xoxo 



oo 



Using (3.13) one can now eliminate D xy in favor of D xy : = (M l ) xy . This produces 
an expression of the same form just with D xy replacing D xy . Next we eliminate D xy in 
favor of D xy , using (3.10). Anticipating D xoxo — — 1/A (which can be shown to hold 
independent of the gap equation) it follows that the saddle point equations SSi/Su x = 
are equivalent to 

-\D XX = 1 , x ^ x , (3.54) 

which is the non-translation invariant precursor of the gap equation used in section 3.3. 
Inserting this into S 2 gives 

S2IU, H] = -\Y J u x u y [Dl y - X 2 Dl xo DlJ + \ £ H xy \D xy . (3.55) 

x,y^x x,y 



So far translation invariance was not presupposed, i.e. we allowed for a nontrivial in- 
dependence in the extremal configurations uj x , x 7^ Xq. We anticipate now from [4] that 
the only solution of (3.54) giving rise to a translation invariant {n x -n y )\^ =QO = —\D xy 
is constant, i.e. 

uj x = const, x 7^ xq . (3.56) 

As seen in section 4.1 then — XD(x) > 1 fixes the constant to be uo^ = ou_(X,V), 
i.e. the unique solution of the translation invariant gap equation — A-D(O) = 1 obeying 
cj_ G (— 4sin 2 7r/L,0). As before we write D_(x — y) for the translation invariant 
propagator with gap u— Further, in the appendix we prove that 



S 2 [u,H] 



> 0, 



(3.57) 



This means that the saddle point is indeed a local minimum of the action (3.51). We 
conjecture that it is in fact a global minimum of the action (3.51) in the domain V{H). 

We now return to the generating functional W[if]. Substituting (3.55) one arrives at 

exp W[H] = exp { - l - £ H xy [XD_(x -y) + l]} 

x,y 

xJ ^J n rfM:r exp { + ^ S u x u y D_ :2 (x,y)^ 

x^x a x,y=£x 

x { 1 + ivTT(-^ + ^) + ((7vTTp)}^ < 3 ' 58) 

Here 

D^ 2 (x, y) := D_(x - yf - \ 2 D_(x - x Q ) 2 D_(y - x ) 2 , (3.59) 



on 



and the terms occurring in the third line are generated by the expansion of 

ex P { " E 7V T l ^-2 Si[u, H] } , (3.60) 

with the half-integer powers of 1/(N+1) omitted. The conversion procedure exemplified 
for S2 is easily seen to generalize to all Si, I > 3. Indeed from (A -1 )^ = —uj Xo + 
(A -1 )'* it follows that in an expansion 

(A 1 ) X0 1 X0 = D x ^ Xo - -j===Q 1 + ^2 ^ N — Yji/2 $ ^xoxo l ' (3.61) 

all coefficients K^ -1 )^ 1 ^]/, I > 2, have the following properties: (i) eliminating D xy in 
favor of D xy = (M~ l ) xy via (3.13), the parameter u drops out. (ii) The result of the 
replacement is an expression of identically the same form, just with D xy occurring for 
each instance of D xy . (iii) In the result one can replace D xy by D_(x — y) and D xy by 
D_(x -y) + \D_{x - x )D-(y - x ). 

The problem of computing VF[if] to order 1/(A^+1) P thus reduces to the algebraic 
problem of computing the Si[u, H] , / < 2p + 2, with the indicated substitutions, and to 
performing Wick contractions with the free propagator 

(u x u y )o = -2A_(x,y) , ^2&_(x,z)D_ :2 (z,y) = 5 X:y , x,y^x . (3.62) 

z 

In fact one has 

A-(x,y) = A_(x - y) , forx,y^x , (3.63) 

where A_ is the translation invariant free -u-propagator defined by J2 Z — z)D-{z — 

On the basis of the general results the if-moments computed with the algorithm de- 
scribed above must be translation invariant and be related to that of the compact model 
by the involution described. As a test of both the algorithm and the involution we now 
compute the two- and the four-point function to sub-leading order. To this end explicit 
expressions for S3 and S4 are needed. 

It is convenient to parameterize them as follows 



+ S^ + = \ [Tr P4 + \Qf + A 2 Q?> + A 3 Qf ] . (3.64) 



on 



Here S z (n) , n = 0,1,2, are the pieces of order n in the source H. The contribution to 



Sj^ coming from [(A 1 )- Eo 1 a;o ]i in (3.61) we denote by Q\ n \ for n = 0,1,2, respectively. 
We need Qf\ Q^p and ■ They come out as 

Qf 1 = ~ ^ u x UyU z D(x-x )D(z-x )D(x,y)D(y,z) , (3.65a) 



g^=2 H xy u z D{x-x )D{z-x )D{y,z) , (3.65b) 

x,y;z^x 

Qi 1 ^- E H xy u z u w D(z-x )D(w,x) 

x,y;z,w^x 

x ^(y-xo)^^, iu) + L>(w-x )5(z, y)] , (3.65c) 

gf= ^ H xy H zw D(x-x )D(w-x )D(y,z) . (3.65d) 

x,y,z,w 

Here we omitted the subscript '— ' on the propagators and we shall continue to do so for 
the remainder of this subsection. 

For Si°\ Si 1] , S£\ S? ] one obtains 

^3°^~g u xU y u z D(x, y)D(y, z)[D(z— x) — 2\D(x— x )D(z — x )] , (3.66a) 

x,y,z^x 

S 3 )= -^ E ^«^(x^)[£)(z- 2 /)-A£)(z-xo)£>(s/-xo)], (3.66b) 



S 1 ! 1 ^- ^ H xy u z u w D(x, w){B(z,w)[D(z-y) - 3XD(z-x )D(y-x )] 

+ D(z,y)[D(w-z) - \D(w-x )D(z-x )}} , (3.66c) 
S?=-^ J2 H xy H zw D(y,z)[D(x-w) - \D(x-x )D(w-x )} . (3.66d) 

x,y,z,w 

The desired correlation functions can now be computed from 

(-5? ) + ^[5? ) ] 2 >o = y E ^. 1 ^,^2 (0) (^i,i/i;^ ) 2/2) ) 

zi, 2/1,^2,2/2 
x,y 



For W 2 ' the xo-dependent terms cancel out algebraically and one finds 

W^\x u y l ;x 2 ,y 2 ) = D(x 1 -x 2 )D(y 1 -y 2 ) + D(x l -y 2 )D(y 1 ~x 2 ) 

-2^D(x l -z)D(y 1 -z)A(z-w)D(w-x 2 )D(w-y 2 ) . (3.68) 

z,w 

The expansion of (— S^) and (S^S^) generates a large number of x -dependent 
terms which simplify after using J2 Z A ( x ~ z )D(z — x ) 2 = 5 XXQ . The results are 

(Si 1] ) = \J2 H xyJ2 A ( z ' w)D(z-w)D(x-w)D(y-z) (3.69a) 

xy z,w 

-A 3 ^2 H xy D(x-x )D(y-x ) ^ A(z, w)D(z-w)D(z-x )D(w-x ) ■ 

xy z,w 

(Si 0) Si% = -\ q J2H xy J2D(z-x)D(z-y) (3.69b) 

xy z 

+A 3 H xy D(x-x )D(y-x ) A (*> w)D(z-w)D(z-x )D(w-x ) . 

xy z,w 

Here 

q := A ( z ) S ^D(z-u)D(z-v)A{u-v)D{u-v) . (3.70) 

2 U,V 

Finally 

Wf 1] (x,y) = -2q^D(x-w)D(y-w) 

w 

+2^D(x-z)D(y-w)A(z-w)D(z-w) , (3.71) 

z,w 

which obeys W^\x,x) = 0. Recall that D{x) here is short for D^(x) = D(x)\ U!=ul _. 
The results (3.71) and (3.68) are manifestly translation invariant and related to those 
of the compact model (see [18, 9]) by the involution described. 



4. Conclusions 

On the level of asymptotic expansions interesting and useful correspondences exist be- 
tween compact and noncompact nonlinear sigma-models. Here we established two such 



correspondences: in perturbation theory, where simply flipping the sign of (3 allows to 
move between the compact and noncompact models, and in the large N expansion, 
where one has to make a detour through the finite volume in order to establish a cor- 
respondence that involves not just flipping the sign of A but also picking the physically 
appropriate solution of the gap equation. 

Taking the sign flip prescription beyond perturbation theory leads from the compact 
model not to the noncompact one, but rather to another compact model, the antiferro- 
magnet. 

In the literature [19] a sign flip rule for the large N coefficients in a continuum formulation 
has been proposed. By section 3.4 the rule is incorrect. In [2] the much weaker claim 
has been made that there exists a scheme in which the coefficients of the large N beta 
function in the compact and in the noncompact models are related by the sign flip 
A — > —A. This statement seems to be in agreement with our findings because it is based 
on re-organizing the perturbation expansion into powers of 1/N, but it remains an open 
question whether the analytic continuation from positive to negative A employed in [2] 
to go from the compact to the noncompact model is justified in the continuum. 
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Appendix A: Positivity of S2 



Here we show that 



S 2 [u,H] 



>0, 



(A.l) 



which guarantees that the saddle point on which the large N expansion is based is indeed 
a local minimum of the action (3.51). 

We begin by noting that the source-dependent term separately obeys 



(A.2) 



for all sources H xy < 0, as —\D(x — y) > 1. It therefore suffices to consider S^w, 0]. 
Next we introduce 



n(p) :=E eip ^) 2 = ^E 



V ^ (E k + u;)(E k _ p 
and claim that in terms of it a sufficient condition for S 2 [u, 0] > is 

n(p)<o, p^o. 



(A.3) 



To see this we prepare 



n(o)>£, 



(A.4) 



(A.5) 



using once more 1 < — XD(x). Then we rewrite S 2 [u, 0] in Fourier space 

S 2 [u,0] = -^^u{pf S m u{q) , 



p,q 



s 00 = -n(o)[\/-A 2 n(o)], 



(A.6) 



'Op 



X 2 v v p , v p := U(p)e- ipxo , 

-5 M yn(p) + \ 2 v*v q , p, g ^ . 

Subject to (A.4) and (A.5) the matrix S has the form of a positive rank one perturbation 
of a positive diagonal matrix, which therefore must be itself positive. 

It remains to show (A.4). We rewrite IT(p) in the form 



n(p) 



E v + 2u 



j + J(p) 



p^O, 



j, \ J_ -^fc + E p _k — E p 
[P) ' V^(E k + Lu)(E p _ k + Lu)' 



(A.7) 



O /I 



using the gap equation. Here J{jp) manifestly has a finite thermodynamic limit. Further 
we claim 

J(p) > for all p , (A.8) 

which implies (A. 4). 

The proof of (A.8) is based on a reorganization of the ci-fold sum in (A. 7) such that 
each term is positive on account of the Lemma below. We begin by noting that J(p) is 
completely symmetric in all arguments p — (pi, . . . ,Pd), and that 

J((27r-pi,pi)) = J(p), p M := (A.9) 

Since J(0) > trivially, it suffices to restrict attention to momenta p ^ with < p\ < 
■ ■ ■ < Pd < 7T- Further, from 

E p . k + £ fe - £ p = 8 cos T Sin ^ Sin ' ( A - 10 ) 



one sees that for L even 

J((7r,...,7r)) = 0, J((0,7r,...,7r))>0. (A.ll) 
Next we insert (A. 10) into (A. 7) to obtain 

~ I^T 2 cos f sin % sin 

" I Z, [sin 2 K _ sin2 f + x,] [sin 2 ^ - sin 2 f + ig ' ( } 

where we set 

X,:=$>n 2 |, F,:=^sin 2 ^^, 2sin|: =v ^. (A.13) 

The result (A.8) now follows from the 
Lemma: For all L > 2 and q e [0,n], q e ?f Z, 

1 2cosf sin^sin(^ - f) 

S(q, «, X, yr-l2^ ^2 ™ _ gin 2 a + X ][ sin 2 ( ™ _ |) _ sin 2 a + y] ^ > ( A ' M ) 



for parameters x, y > and < a < 7r/L. 



O IT 



In the application to (A. 12) only the bound on a needs to be verified. From (3.28) we 

have however 

a / 4 7r 

sin— < 4/ — sin — . (A. 15) 

2 \2d+l 2L K ' 

Proof of the Lemma. Trivially S(0, a, x, y) > and S(tt, a, x, y) — for L even. It thus 
suffices to consider q e (0, n). We make use of the following summation formulas: 



1 1 _ 2coth^ 

sin 2 ^ + sh 2 § ~ sha 

n=0 2 2 



a^O, (A.16) 



and for g e ^ Z and a ^ o ^ 0: 



L-l 



L ^ [sin 2 f + sh 2 f][sin 2 (^ - f ) + sh 2 |] (A ' 1?) 



2 



sh^sh^shash6 



shf(a + b) sh(a + 6) sh|(a - 6) sh(a - b) 



+ 



ch(a + 6) — cos q ch(a — 6) — cos g 



For a = b ^ the result of the summation is obtained from the rhs of (A. 17) by omitting 
the second term in square brackets and specializing the remainder to a = b (i.e. no pole 
occurs for a = b and q — 0). Both (A.16) and (A. 17) can be obtained from 



lEl^^. (A.18) 
L n=0 e l - a 1 - a 



In the application to S we set sh | = x — sin 2 |, sh | = y — sin 2 |, and pick the roots 
a, 6 > for x, y > sin 2 |. In a first step one obtains (also using the Id version of (A. 10) 
backwards) 



S = 



2coth^ 2cothf 2(sh 2 f + sh 2 | + sin 2 f) 



sha 



+ 



shfe 



sh^fshf shashfr 



x 



shf(a + 6) sh(a + 6) shf (a - b) sh(a - 6) 



ch(a + 6) — cosg 
This can be rewritten in the form 



ch(a — 6) — cos q 



(A.19) 



cos 2 f 



chf chf sh^sh^ 



shf (a + b) sh^ shf (a - 6) sh^ - 



sh 2 ^ + sin 2 f 



sh 2 ^ + sin 2 § 



(A.20) 



Since the function 



t(z) := 



sh^f shf 



sh^ 



sin 



2 <? 



(A.21) 



is positive and monotonically increasing on IR + it follows that for a, b > and q G (0, ir) 
the rhs of (A. 20) is nonnegative. 

Next consider the case where < 8 := ib < n/L and a > 0. The expression (A. 20) can 
in this case be written as 



4 cos 2 f 
S = 2 



coth^f si - cot^s 2 



ch| cos | (cha cos (3 — cos q) 2 + sh 2 a sin 2 (3 

i ° Pi i ^ l 

Si = sh- cos — [cha — cosp + 1 — cosgj , 
s 2 = ch- sin ^[cha — cos/? — 1 + cosg] . 



(A.22) 



Further 



La L(3 . 1 n . , a /3 

coth— si — cot — s 2 = (cha — cos p) ch- cos — 



2 J 
+ (1 — cosg) 



La a L/? /3 

coth— th- — cot — — tan 



La a (3 L8 . 3 a 

coth — sn- cos — h cot — sin — cn- 
2 2 2 2 2 2J 



2. 

(A.23) 



The second term in (A.23) is manifestly positive. For the first term we note that 8 i— > 
cot tan | maps [0, n/L] bijectively onto [1/L, 0], while a h- > coth^th| maps [0, oo) 
bijectively onto [1/L, oo). Hence 



La a L8 (3 

coth — th cot — tan — > , 

2 2 2 2 ~ 



(A.24) 



which shows that S in (A.22) is nonnegative. 

It remains to consider the case < 8 := ib < n/ L, and < a := ia < ir/L, (which may 
differ from the a used in the formulation of the Lemma). Specialization of (A. 20) gives 



cos 2 | 



cos f cos f sin ±f sin M 



sinf(a + /3) sin ^ sin f (a - 8) sin 2=2 



sin 



2 9 



sin 



2 a+P 



sin 2 | 



sin 



2 Q-/3 



(A.25) 



Since x \— > —t(ix) (with t from (A. 21)) is again positive and strictly increasing for 
x G [0, 7r/L] and g G (0, 7r), it follows that the rhs of (A.25) is non-negative for the range 
of a,j3 considered. This concludes the proof of the Lemma. 
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